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We investigate the phase diagram of the so-called Polyakov-Nambu- 
Jona-Lasinio model at finite temperature and nonzero chemical potential 
with three quark flavors. Chiral and deconfinement phase transitions are 
discussed, and the relevant order-like parameters are analyzed. A special 
attention is payed to the critical end point (CEP): the influence of the 
strangeness on the location of the CEP is studied; also the strength of the 
flavor-mixing interaction alters the CEP location, once when it becomes 
weaker the CEP moves to low temperatures and can even disappear. 
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Understanding the QCD phase structure is one of the most important 
topics in the physics of strong interactions. The developed effort on both 
the theoretical point of view (by using effective models and lattice calcula- 
tions) and the experimental point of view (it is one of the main goals of the 
heavy ion collisions program [1]) has proved very fruitful, shedding light on 
properties of matter at finite temperatures and densities. 

Confinement and chiral symmetry breaking are two of the most impor- 
tant features of quantum chromo dynamics (QCD) which is the fundamental 
theory of strong interactions. Its basic constituents are quarks and gluons 
that are confined in hadronic matter. At high temperatures and densities 
hadronic matter should undergo a phase transition into the quark-gluon 
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plasma (QGP). A challenge of theoretical studies based on QCD is to pre- 
dict the equation of state, the critical end point and the nature of the phase 
transition. 

The fundamental theory of QCD is known to be highly difficult for the 
analytical analysis because of the strong coupling regime. Meanwhile, as 
this regime triggers the most interesting phenomena in QCD, it has become 
a common practice to replace the low-energy QCD by some effective mod- 
els. Examples of such descriptions include NJL type models which have 
been developed providing guidance and information relevant to observable 
experimental signs of deconfinement and QGP features. 

NJL type models, that take into account only quark degrees of freedom, 
give the correct chiral properties; static gluonic degrees of freedom are then 
introduced in the NJL Lagrangian through an effective gluon potential in 
terms of Polyakov loop with the aim of including features of both chiral 
symmetry breaking and deconfinement. 

Our calculations are performed in the framework of an extended SU(3) j 
PNJL Lagrangian, which includes the 't Hooft instanton induced interaction 
term that breaks the Ua(1) symmetry, and the quarks are coupled to the 
(spatially constant) temporal background gauge field $ [21 [3]: 

1 8 

C = q(i-fDp-m)q+-g S £ [ ( q X a q f + ( qi 75 A a q f ] 

Z a=0 

+ g D {det[q(l + j 5 )q]+det[q(l- l5 )q]}-U($[A},$[A];T). (1) 

The covariant derivative is defined as = — iA^, with A^ = 5qAq 
(Polyakov gauge); in Euclidean notation Aq = —1A4. The strong coupling 
constant g is absorbed in the definition of A >Jj (x) = gA£(x)^f- , where A% is 
the (SU(3) C ) gauge field and A a are the (color) Gell-Mann matrices. 

The effective potential for the (complex) field adopted in our parame- 
trization of the PNJL model reads: 



W($,$;T) _ a{T) 
T 2 ~ 2~ 



$$ + 6(T)ln 1 - 6$$ + 4($ 3 + $ 3 ) - 3($$) 2 ,(2) 



where 



a (T) =a + a 1 + a 2 ^) and 6 ( T ) = 6 3 



(3) 



The parameters of the effective potential U are given by ao = 3.51, 
ai = —2.47, 02 = 15.2 and ^3 = —1.75. These parameters have been 
fixed in order to reproduce the lattice data for the expectation value of the 
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Fig. 1. The phase diagram in the SU(3) PNJL model: the location of the CEP is 
found at T CEP = 155.80 MeV and fi CEP = 290.67 MeV (see details in text). 

Polyakov loop and QCD thermodynamics in the pure gauge sector. When 
quarks are added, the parameter To, the critical temperature for the decon- 
finement phase transition (that manifests itself as a breaking of the center 
symmetry) within a pure gauge approach, was fixed to 270 MeV, according 
to lattice findings. This choice ensures an almost exact coincidence between 
chiral crossover and deconfinement at zero chemical potential, as observed 
in lattice calculations. 

The parameters of the NJL sector are: m u = rrid = 5.5 MeV, m s = 140.7 
MeV, g s A 2 = 3.67, g D A 5 = -12.36 and A = 602.3 MeV, which are fixed to 
reproduce the values of the coupling constant of the pion, f n = 92.4 MeV, 
and the masses of the pion, the kaon, the r\ and r)', respectively, = 135 
MeV, M K = 497.7 MeV, M v = 514.8 MeV and M v > = 960.8 MeV. 

The inclusion of the Polyakov loop effective potential Z//(<3?, T), that 
can be seen as an effective pressure term mimicking the gluonic degrees of 
freedom of QCD, is required to get the correct limit. Indeed in the NJL 
model the ideal gas limit is far to be reached due to the lack of gluonic 
degrees of freedom. 

The location and even the existence of the CEP in the phase diagram is a 
matter of debate. While different lattice calculations predict the existence of 
a CEP [4j, the absence of the CEP in the phase diagram is still a possibility 
as was seen in some lattice QCD results [5], where the first order phase 
transition region near [i = shrinks in the quark mass in fi space when \jl 
is increased [5]. 

In Figure [T] is presented the phase diagram in the PNJL model. As the 



4 



sqm2011'PedroCosta printed on December 30, 2011 




H (MeV) H (MeV) 

Fig. 2. Left panel: the phase diagram in the SU(3) PNJL model. The solid lines 
represent the first order phase transition, the dotted line the second order phase 
transition, and the dashed line the crossover transition. Right panel: the phase 
diagram and the "line" of TCPs for different values of m s (the dotted lines are 
just drawn to guide the eye); the TCPs in both figures are obtained in the limit 
m u — rrid = and m s ^ 0. 



temperature increases the chiral transition is a first order one and persists 
up to the CEP. At the CEP the chiral transition becomes a second order 
one. The location of the CEP is found at T CEP = 155.80 MeV and \£ EP = 
290.67 MeV [p% EP = 1.87po)- For temperatures above the CEP there is a 
crossover whose location is calculated making use of d 2 (qq) jdT 2 = 0, i.e. 
the inflection point of the quark condensate (qq). 

The transition to the deconfinement is given by d 2 &/dT 2 = 0, and is 
represented by the magenta line. The surrounding shaded area that limits 
the region where the crossover takes place is determined as the inflection 
point of the of the susceptibility d&/dT. 

Due to the importance of the location of the CEP from the experimental 
point of view, let us investigate the influence of several parameters which 
can lead to a significant change in the CEP's localization. We study the 
influence of strangeness on the location of the CEP (or tricritical point 
(TCP) in the chiral limit). 

A first point to be noticed is that in the PNJL model, when the full 
chiral limit is considered (m u = rrid = m s = 0), the phase diagram does not 
exhibit a TCP: chiral symmetry is restored via a first order transition for 
all baryonic chemical potentials and temperatures (see left panel of Figure 
[2|). On the contrary, in the chiral limit only for the SU(2) (light) sector 
(m u = rrid = 0, m s ^ 0) a TCP can be found [6]. Both situations are in 
agreement with what is expected: the chiral phase transition at the chiral 
limit is of second order for Nf = 2 and first order for Nf > 3 [7j . 
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To study the influence of strangeness on the location of the critical 
points, we vary the current quark mass m s , keeping the SU(2) sector in 
the chiral limit and the other model parameters fixed. In Figure [2] (left 
panel) we plot the phase diagram of the model in the (T, fj,) plane, for 
various values of the the current quark mass m s . 

The pattern of chiral symmetry restoration via first order phase tran- 
sition remains for m u = = and m s < mf lt . The value for mf lt is a 
subject of debate; we found mf lt ~ 9 MeV in our model, lower than lattice 
values [8] and half of the value obtained in NJL model (mf" lt = 18.3 MeV 
[6]). When m s > mf* 1 , at \i = 0, the transition is a second order one and, 
as /i increases, the line of the second order phase transition will end in a 
first order line at the TCP. Several TCPs are plotted for different values 
of m s in the right panel of Figure [21 As m s increases, the value of T for 
this "line" of TCP's decreases as fi increases getting closer to the CEP and, 
when m s = 140.7 MeV, it starts to move away from the CEP. The TCP for 
m s = 140.7 MeV is the closest to the CEP and is located at ji TCP = 206.95 
MeV and T TCP = 187.83 MeV. If we choose m u = ^ 0, instead of a 
second order transition we have a smooth crossover for all the values of m s 
and the "line" of TCPs becomes a "line" of CEP's. 

Also the change of the Ua(1) anomaly strength has a strong influence 
on the localization of the CEP in the (T, fx) plane. In Figure El we show the 
location of the CEP for several values of go compared to the results for go Q , 
the value used for the vacuum. As already pointed out by K. Fukushima in 
[2], we also observe that the location of the CEP depends on the value of go- 
In fact, our results show that the existence or not of the CEP is determined 
by the strength of the anomaly coupling, the CEP getting closer to the \x 
axis as go decreases. 

As a conclusion, we investigated the phase diagram of the so-called PNJL 
model at finite temperature and nonzero chemical potential with three quark 
flavors. Chiral and deconfinement phase transitions are discussed, and the 
relevant order-like parameters are analyzed. 

The chiral phase transition is a first order one in the chiral limit. Work- 
ing out of the chiral limit, at which both chiral and center symmetries are 
explicitly broken, a CEP which separates first and crossover lines is found, 
and the corresponding order parameters are analyzed. 

A special attention is payed to the critical end point. We studied the 
influence of strangeness on the location of the critical points varying the 
current quark mass m s , keeping the SU(2) sector in the chiral limit and the 
other model parameters fixed. When m s > mf , at /i = 0, the transition is 
of the second order and, as \x increases, the line of the second order phase 
transition will end in a first order line at the TCP. If we choose m u = ^ 0, 
instead of a second order transition we have a smooth crossover for all the 
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Fig. 3. Dependence of the location of the CEP on the strength of the 't Hooft 
coupling constant go- 

values of m s and the TCP becomes a CEP which depend strongly on the 
the value of m s . 

We also analyzed the effect of the anomalous coupling strength on the 
location of the CEP. We showed that the location of the CEP depends on the 
value of go and, as this strength of the flavor-mixing interaction becomes 
weaker, the CEP moves to low temperatures and can even disappear. 

Work supported by Centro de Ffsica Computacional and F.C.T. under 
Project No. CERN/FP/116356/2010. 
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